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We report an experimental, numerical and theoretical study of the motion of a ball on a rough 
inclined surface. The control parameters are D, the diameter of the ball, 9, the inclination angle of 
the rough surface and Eki, the initial kinetic energy. When the angle of inclination is larger than 
some critical value, 9 > 9t, the ball moves at a constant average velocity which is independent of 
the initial conditions. For an angle 9 < 9t, the balls are trapped after moving a certain distance. 
The dependence of the travelled distances on Eki, D and 9. is analysed. The existence of two kinds 
of mechanisms of dissipation is thus brought to light. We find that for high initial velocities the 
friction force is constant. As the velocity decreases below a certain threshold the friction becomes 
viscous. 



I. INTRODUCTION 



Essentially due to the non linear and dissipative character of their interactions, flowing grains exhibit a range of 
complex and fascinating behaviour such as density waves avalanches arching and segregation Q. The 
energy exchange among the grains and between grains and walls are particularly important, but not well understood. 
The solid friction force itself is still an open problem dating back to the work of Coulomb. In addition, the collective 
effects of this solid friction in a quasistatic deformation of a granular solid are not well understood ^ . To take into 
account the energy loss during collisions, we use only a very empirical parameter -the restitution coefficient- even if 
experiments Q] have shown that this coefficient is more complex than expected. 

We have performed several experiments |^ on the motion of balls on rough inclined planes. One of our most 
important results is that in a regime characterized by an average constant velocity V , the macroscopic friction force is 
"viscous" , i.e. proportional to V . In the D — 9 parameter space, where D is the diameter of the moving ball and 9 the 
inclination angle, the constant velocity regime, which we call regime B, is sandwitched between a decelerated/trapping 
regime A (for small D and/or 9) and a regime C in which the ball jumps on the plane, and eventually can reach a 
chaotic motion In this paper, we focus on how the ball losses its energy in the decelerated/trapping regime, A. 

After a short description of the experimental procedure, we present results of the stopping distance of balls reached 
with low initial kinetic energy, as well as results from an approximate two dimensional stochastic model. We discuss 
energy dissipation in section |I|, and in the last section we propose a phenomcnological model and compare it with 
numerical simulation. 



II. EXPERIMENTAL PROCEDURE 



The experimental setup used in this work has been described in detail elsewhere 11 . A rough plane is constructed 
by sticking particles (see below) on an adhesive surface. This surface is placed on a thick glass plane supported by 
a metal frame to prevent warping and a jack that allows the accurate adjustment and measurement of the angle of 
inclination, 9, which is one of the control parameters in our experiment. The rough surface is made of a monolayer of 
sifted rolled sand grains with mean radius, r, 0.2mm < r < 0.25mm. These grains are glued to the adhesive surface 
in such a way as to obtain a homogeneous geometrically disordered rough surface with a surface fraction compacity 
close to 0.8. We used this surface to study the motion of steel spheres of diameter 1.6mm < D < 10. 3mm. 

In order to control the initial kinetic energy supplied to the ball (i.e. the initial speed with which the ball hits the 
rough surface) a thin smooth plastic sheet is placed on the rough surface in the region where the ball is launched. A 
set of parallel and equally spaced lines is drawn on the sheet perpendicular to the direction of maximum slope of the 
plane. The first line is drawn at a distance equal to the smallest rolling ball radius, from the edge of the sheet nearest 
to the rough surface. To study the behavior of rolling balls with very low initial kinetic energy, Eki ~ 0, we release 
them from this first line. To increase the initial Eki, we simply release the ball from one of the other lines traced 
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higher up on the plane. It is then trivial to see that Eki oc Xi, where Xi is the distance travelled over the smooth 
surface. 

The experimental procedure is as follows. The inclination angle of the plane is fixed at a desired value, and a 
straight barrier is placed at the location which gives the desired initial kinetic energy, Eki- Twenty balls of a given 
radius, i?, are positioned on the barrier on the uphill side and not touching each other. Removing the barrier, the 
balls are accelerated as they move the same distance, Xi, on the smooth surface, arriving at the rough surface with 
the same initial kinetic energy, Eki- Since we have set the inclination angle, 9, in the range corresponding to the 
pinning regime. A, the released balls will come to a stop on the rough plane at various distances Li. The susbscript i 
denotes a particular pinned ball. This procedure is repeated 35 times with the same experimental settings to collect 
statistics. We then repeat the experiment for different values of Eki, R and 9, always keeping 9 in the pinning regime. 
The uncertainty in length measurements is of the order of ±1 mm and the number of balls for each set of parameters 
values investigated is large enough to have good statistics. 



III. STOPPING DISTANCE OF BALLS WITH LOW INITIAL KINETIC ENERGY 

With the method described above, we measured the distance, Li, travelled on the rough surface by each ball and 
for a wide range of control parameter values (1.6mm < D < 10.6mm and 1.5° < 6* < 12.5°). In all cases Xi is taken 
as small as possible, around 2 mm. In Figs. ^ (a, b) typical Li distributions are shown for D — 2mm and 9 — 3° and 
5° respectively. It can be seen that these distributions peak for a travelled distance, Li, of the order of the diameter 
of the moving ball (2mm). The distributions decay for larger values of the travelled distance, the decay being faster 
for the smaller 9. 

We found that the decay of these distributions for the larger values of Li is well fitted by an exponential 

N = Noe-°'^\ (3.1) 

The solid lines in Figs. |l|(a) and (b) are fits obtained with a — 2.9cm^^ and 0.5cm^^ respectively. Similar exponential 
distributions were found for the stopping distance by Riguidel et al, but working under experimental conditions 
corresponding to regime B, where the moving ball reaches a steady, albeit fluctuating, velocity. In this regime, the 
balls are occasionally stopped by large holes on the surface. Their results were well fitted with a varying with 9 as: 

a = exp {-aD^ sin^ 9) , (3.2) 

where a is a constant. 

In order to complete and verify this law, we have performed a systematic study of the variations of a with 9. In 
Fig. ^ we show a as a function of 9 for D = 2mm. Similar results were found for all the D values studied. Two 
different dependences are observed: A linear variation with large negative slope for inclination angles lower than 4° 
followed by a crossover to a dif fere nt smoothly decreasing curve tending to zero for angles larger than 5°. Good 



agreement was found using eqn. 3.2 for this last part of the curve, as shown in the figure (solid line). However, it is 
clear that eqn. 3.2 is not valid for 9 values smaller than 4°. This clearly demonstrates the existence of a change in the 
trapping mechanism of the rolling ball when the inclination angle goes from a value lower than 4° to a larger one. To 
find the angle at which the first trapping mechanism disappears, we look for the intercept of the straight line in Fig. |^ 
with the x-axis. We define this transition angle, 9t, a = as the angle for which the trapping mechanism of regime 
A has disappeared. In Fig. |^, 9t = 4.6°. This value is very close to that found previously for the angle 9ab 
at which the transition between the decelerated regime (regime A) and the mean constant velocity one (regime B) 
occurs (for the same D). In the same way, we can find the angles 9t corresponding to a = 0, for all studied D values. 
These 9t angles can be compared with the 9ab values previously reported ||l^ . Very good agreement is obtained for 
large D values, and slight differences for smaller ones. 

We can therefore conclude that there are two trapping mechanisms: One due to large wells (disorder) in the surface, 
and the other one to dissipation (controlled to a large degree by the relative smoothness seen by the rolling ball). 
The measurement of the dependence of a on described above, provides a more physical and precise criterion for 
determining experimentally the transition line between these two regimes. 

The fact that in Fig. |^ the transition between the two observed behaviours according to 9 values occurs smoothly 
indicates that near the transition, both trapping mechanisms are important. 

In fact, this experimental distribution of the stopping distances and their exponential fits was predicted using a 
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simple two dimensional^ stochastic model. The above experimental verification came later. In this model one starts 
with Newton's equations of motion for the moving balls which are then simplified by making some approximations 
motivated by the geometry of the collisions and the properties of the dissipation. The resulting model describes 
well the statistical properties (i.e. averages and distributions) of physical quantities of the moving balls. For details 
of the model see reference ||l2|. In Fig. |^ we show the a versus 9 plot obtained from the stochastic model. We 
see that it has the same form and functional dependence as that found in the experiments. The actual values of a 
and 9 given by the model differ from the experimental ones. The model is too simple to give accurate quantitative 
predictions, for example rotation is ignored. However, the excellent qualitative agreement between the predictions 
and the experimental results means that the role of geometry, which is emphasized in the model is indeed crucial. 

We define the median, L*, of the stopping distance distribution as the length for which 50% of the balls are trapped 
on the rough surface. Since we found exponential distributions in all cases we have studied, we can write L* as: 

L* = (3.3) 

a 

As seen in Fig. § a = — asin(6') + 6, and from the definition of 9t, aifix) — 0, we see that ae^Or — a{9 — 9T) cos{9t), 
and therefore Lg^g^ cx j^j^g^ j . Thus, L* diverges for ^ 0^ in a way reminiscent of a phase transition with a critical 
exponent equal to unity. 

In order to verify this behaviour we measure L* for the distributions corresponding to different pairs of control 
parameter values {9,D). In Fig. ^ the dependence of L* on 9 is shown for three D values. For D = 10.3mm the 
transition between regimes A and B occurs at 6* = 2.5° and for D — 3mm at = 4.4°. In both cases the very rapid 
increase of L* as — > 9t from below, demonstrates an approximate divergence consistent with the above discussion. 
Let us recall, that as — > 0t we have two competing pinning mechanisms and therefore no true divergence. This 
"divergence" was based on the idealization that as we approach 9t from below, pinning in regime A completely 
disappears. For the smallest D value shown in the figure {D = 1.6mm), the agreement between the proposed 
functional dependence of L* on 9 and the experimental results is also clear, even if the A-B transition is not yet 
reached. 

IV. ENERGY DISSIPATION IN REGIME A 

In this section we discuss energy dissipation in regime A, i.e. small angles of inclination, 0, where the balls always 
come to rest. With this goal in mind, balls were released from the various lines marked on the plastic sheet, as 
described in section ||. In this way, controlled initial kinetic energies were supplied to the balls, and their stopping 
distances were studied. 

We will first show that the transition A-B is not affected by the initial kinetic energy. This is in agreement with 
previous characterization of regime B as that interval of inclination angles for which the balls reach an average 
steady state velocity independent of the initial velocity (or kinetic energy). |p2|-p^ We then analyze, for a very small 
inclination angle {9 = 2°), the stopping distance distributions for different initial kinetic energies and several D values. 
After that, we investigate the dependence of the mean distance travelled by the ball with its initial velocity. Finally, 
we propose a phenomenological interpretation of the experimental results obtained. 

A. Experimental results 

To study the influence of the initial velocity {Vi) on the transition between regimes A and B, 100 balls of diameter 
D = 3mm were released from Xi — 0.3 cm {Vi ~ 5 cm/s) and then from Xi = 20 cm (Vi ~ 45 cm/s). The stopping 
distance of each ball was measured, and L* calculated. This procedure was repeated for ten different inclination 
angles of the rough surface, 2.7° < 6* < 6°, in order to insure the change of the dynamical regime. 

We have shown before that, with its rapid increase ("divergence"), L* itself characterizes the transition. In Fig. ^ 
we show this divergence of L* as a function of the inclination angle for the two values of the initial velocity, Vi = 5,45 
cm/s. It is clear that the A-B transition occurs at the same angle for both curves independently of Vi. This is so even 
though the initial velocity has a great influence on the travelling distance in regime A. 



^The experiments are, of course, three dimensional: The two directions defining the rough plane, and the direction perpen- 
dicular to it. 
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We now discuss the evolution of the distribution of stopping distances as a function of the initial kinetic energy of the 
balls. The experimental procedure is as follows: The control parameters (9 and D) were fixed at values corresponding 
to regime A and a large enough number of balls were released at different Xi. The distance travelled by every ball 
before being trapped was measured. We show the distributions for the stopping distances in Fig. || for D = 10.3mm, 
= 2° and four values of Vi. The distributions are essentially the same, except that the location of their centers 
moves towards increasing values of the travelled distance as the initial velocity is increased. The dispersion is very 
small. 

We observe experimentally that those balls which begin their motion with larger initial velocity travel a longer 
distance before getting trapped. Moreover, balls released from the same Xi (i.e. the same initial velocity) move 
through almost the same distance on the rough surface before getting trapped over a distance of a few centimeters: 
The dispersion of the distribution of stopping distances is very small. This leads us to postulate that the ball does 
not get trapped unless its velocity is under some threshold value. To arrive at this value, it must first travel a certain 
distance on the rough surface to dissipate enough energy. Using a video camera and image processing, we were able 
to evaluate this threshold value. We launch several balls {D = 6mm, 9 = 2°) with different initial kinetic energies 
(0.5cm < Xi < 18.5cm) and for each ball, the velocity is measured every acquired frame (15 frames by second). Balls 
travel a distance of about 15 cm before getting trapped. Our measurements indicate that balls do not get trapped 
while their speed exceeds 3 cm/s. Note that the velocity below which the ball can be trapped is independant of the 
initial velocity^. Since the distributions of stopping distances are rather narrow and symmetric, the average value L 
is a good variable to use to characterize the energy loss. We show in Fig. ^ a plot of X as a function of Xi, for 9 = 2° 
and D = 10.3, and 4.7mm. 

The first thing to notice is that L increases with the initial kinetic energy (oc Xi) of the balls. We also see that for 
D = 10.3mm and Xi > 10 cm, this dependence is linear. On the other hand, for 0cm < Xi < 10cm, L increases with 
Xi following another law. Note that for Xi = Q, L must be 0. The same behaviour is seen for D = 4.7mni but, in 
this case, the crossover between the two dependences is around Xi = 2 cm. 

The presence of two different behaviours for the stopping distance as a function of the initial kinetic energy suggests 
a change of the nature of the friction force between the two regimes. At the moment we cannot characterize these 
forces any further. However, in previous work ]^ Jl^ , p^ , |lT[| , which concentrated on regime B, friction mechanisms 
were extensively studied experimentally, numerically and theoretically. Only three types of friction force are possible: 
F = Ki, F = K2V and F = K^V^ . In the following section we use these experimental results to propose a 
phenomenological model with few parameters which reproduce the experimental results. 

B. Physical model 

Our objective is to describe the dynamics of a ball moving down a rough slightly inclined surface. Taking into 
account the experimental results just described, we can assume that two different types of friction forces exist. Which 
type of force enters into play depends, among other things, on the initial velocity of the ball. It is clear that when 
moving down the plane, a ball loses its initial energy by collisions and friction with the surface grains and finally is 
trapped. However, not much is known about the way in which this occurs. In particular, nothing is known about the 
average velocity as a function of the distance travelled. So, the main assumptions of our model are: (a) it is possible 
to define such a velocity function, V{D), for all balls, and (b) to get trapped, V must be smaller than a certain 
threshold value Vmin- In other words, we are assuming that to decrease its velocity from a value Vb to a smaller one 
Vi, a ball has to travel the same distance Aqi independent of its initial velocity. 

As for V{D) we know that: V{0) = Vi and we assume that V{L) < Vmin where Vi is the initial velocity of the 
ball and L is the corresponding trapping distance. Clearly, since the physical quantities of interest are fluctuating 
quantities given by some distribution, the quantities entering in the model are statistical averages. 

With the above considerations, we can write: 

Lo = Li+Xqi, (4.1) 

where Lq and Li are the average stopping distances of balls released with initial velocities Vq and Vi < Vo, respectively, 
and Aoi is the mean distance travelled by the first set of balls as their velocity decreases from Vq to Fi . Differentiating 
this equation with respect to V'^, we obtain: 



■^The same experiment was performed for D = 6mm, and 6 = 2.85° the threshold obtained in this case is 6.5 cm/s. 
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dX{V) _ dL{V) 



(4.2) 



The right hand side of eqn. 4.2 may be evaluated from the experimental results presented in the previous section 
(Fig. As we have already mentioned, for large enough initial kinetic energies, the dependence of L on V'^ is linear. 

Therefore, eqn. 4.2 leads to: dV^ /dX{V) — constant. This is the energy gradient and in consequence must be equal 
to the force. We thus find that, at large enough velocities, a ball moving down a rough surface suffers a constant 
friction force, F — mK , and that the experimental determination of L{Eki) provides a way to evaluate it. 

On the other hand, for smaller initial kinetic energies we have found a different relation between L and which, 
in fact, implies another mechanism for energy dissipation. In this case, we assume a viscous type friction force (i.e. 
proportional to the velocity). For each pair of values {9, D), we define a velocity, VJ, for which the crossover between 
the two types of frictional forces occurs. 

In other words, we propose that a ball of diameter D, moving down a rough surface made of grains of mean radius 
r, and inclined an angle 9, suffers a constant dissipation force {F — mK), while its velocity V is larger than Vi{9, D); 
a viscous type dissipation force {F — niaV) when its velocity is between Vi{9,D) and Vmin{9, D), and that it gets 
trapped on the rough surface only if V < Vmin {9, D) . 

To verify these assumptions, we calculate the distances travelled given by the above model and compare them to 
the experimental ones. If a ball is released with a velocity V > Vi, the distance it travels according to this model is 
easily calculated to be 



On the other hand, if the release velocity, V, is less than Vi, the distance travelled is 

.-07)4(1-^ 



(4.3) 



(4.4) 



We arrive then, at two equations with four parameters: a,K (which characterize the friction forces), Vmin and Vi. 
Us ing eqn. 4.2, K{D) was calculated by fitting to the linear part of the L versus curves, as explained just after 



Eqn. 4.2. Since the slopes of the linear regions in Fig. ^ depend on _D, K is then a function of D. To determine a 
and Vmin we fit Eqn. 4.4 to the curved regions (small Xi) of Fig. ^. However, in this region, we have enough data for 
such a fit only for D — 10.3mm, which is what we use. We note that for small Xi (small initial velocity) the stopping 
distance appears to be roughly the same for all D values we studied. We, therefore, take the values determined for 
D = 10.3mm as constant for all D. 

Finally, the values of Vi for the different D values were evaluated as follows. The constant (i.e. independent of V) 
term of eqn. |4.3|, 



•' 2K a \ Vi 



(4.5) 



can be obtained directly, for each D, from the intersections of the fitted straight lines of L versus V^ with the 
vert ical axis. The solution of eqn. 4.5 gives two values of Vi{D), which are the two intersection points of equation 4.3 



and 4.4, The smaller one is taken as the physical one corresponding to a passage from a more dissipative way to a 
less one. Fig. || shows the result of this analysis for all the studied values of D. Here we plot the experimental average 
stopping distan ce, L , as a fu nction of the initial velocity (not Xi oc V^). The solid lines show the calculated values 
using equations 4.4 and 4.3, It can be seen that the agreement is very good. Note that since a and Vmin have been 
taken as constants for all D, we only have two parameters left for adjustment to fit the data for the various D values. 
So, the very good agreement with experiments is obtained by tuning only two parameters. 



I>(mm) 


2 


3 


4.7 


6.3 


7.1 


10.3 


K{m.s-'^) 


4.9 


2.6 


1.2 


0.8 


0.7 


0.5 


Vi{cm.s-') 


1.2 


2.8 


7.3 


9.1 


10.8 


20.7 



Table I. Values of the different parameters, for 9 = 2° 
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The values of the parameters are given in table The value obtained for Vmin, for the largest D, is consistent with 
zero: it is not possible to give a L value acceptable when it is smaller than D, due to the measurements precision 
and the width of the distributions. We have obtained a = 0.03s~^ for the prefactor in the viscous force while K and 
Vi values range between 500cm. and 50cm. and 1 cm/s and 21 cm/s, respectively, for values of D increasing 
from 4 to 20.6. It can be seen that the maximum value of the viscous deceleration, corresponding to the largest Vi, 
is around 0.6cm. s^^ . 

This means that a ball with D — 10.3mm which begins its motion over the rough surface with, for example, a 
velocity of 50 cm/s, "feels" a constant deceleration of 50cto.s~^ until it reaches a velocity of 21 cm/s. At that 
moment, the force abruptly becomes viscous and the deceleration decreases to 0.6cto.s~^. The velocity then continues 
to decrease and so does the friction force. When the velocity reaches the minimum value, Vmin, the ball gets trapped. 

This behaviour is easily motivated physically. When the ball is launched with a high velocity, collisions with the 
bumps on the surface cause it to undergo rather large bounces. The time of flight, dt, of these (ballistic) bounces is 
determined primarily by V± (the velocity normal to the surface), dt oc Vj_, and therefore the frequency of collisions is 
proportional to V£^. With each collision the ball losses an amount of energy proportional to V±, due to the coefficient 
of restitution. Therefore the energy lost per second, the product of these two quantities, is constant, i.e. a constant 
friction force. More elaborate calculations find similar results. |^ 

At smaller velocities, the bounces are not high, the ball probes the geometry of the surface and the motion is a 
mixture of bouncing and rolling. The time of flight is, therefore, more complicated. The motion closely resembles 
that in the constant velocity regime thus giving viscous friction. The difference between this motion and the constant 
velocity regime, is that here the energy gained by the ball moving down the slope cannot compensate for the energy 
lost. The ball eventually gets pinned. 

Vi appears like a limit velocity at which the ball can just "fly-over" the grains constituting the surface. 



V. NUMERICAL SIMULATIONS 



To investigate in more detail our conjectures concerning the motion of the particle we performed numerical simu- 
lations of the system. The motion of the particle was simulated using soft sphere molecular dynamics (for details see 

f). The sphere moves on a plane configuration scanned in from one of those used by Riguidel in his experiments 
As material parameters we use a normal coefficient of restitution e„ = 0.6 and a coefficient of friction oi ^ = 0.13. 
First we checked that we find the same global behaviour as in the experiments. Figure ^ shows the mean stopping 
distance, as averaged over the launching of 60 balls per starting velocity as a function of the initial kinetic energy. We 
find the same behaviour as in experiments - a linear region for higher starting velocities, and a not so well defined 
different shape of the curve for smaller starting velocities. 

We now want to check whether the ball really covers the largest part of L in large jumps. Figure |l^ shows the 
distance covered between jumps as a function of time for different starting velocities. Obviously, the first part of the 
motion consists of very wide jumps (covering a few particle diameters at high starting velocities), but then drops 
very rapidly to much smaller distances. The times between collisions suddenly exhibit the very regular behaviour 
also observed in simulations of the steady state motion |l^,|l^: a number of small jumps, in the course of which all 
normal velocity with respect to the ball on the plane is lost, followed by a rolling over the rest of the ball. With the 
small angles of inclination and velocities in this case, this rolling starts very early on the ball, as is obvious from the 
long distances between collisions corresponding to them. In Fig. ^ the total distance covered in the same runs as in 
Fig. |l^ is shown. Here, it becomes obvious that indeed, most of the L is due to the large jumps at the beginning of the 
motion, and we see a very clearly defined crossover in the friction force. Only the curve corresponding to Fig. px|(a) 
does not exhibit this crossover. The reason is that the initial jumps were already not much longer than a particle 
diameter, since the initial velocity was quite low. Thus, the discrepancy between the two regimes of the motion is not 
very strong. Besides, the starting velocity corresponding to Fig. |l^(a) is quite close to the lower limit of the linear 
region in Fig. 0, i.e. close to a different type of behaviour. 

From Fig. M it is also obvious that the stopping time is very similar in all cases, but that the onset of rolling 
appears a bit later for higher starting velocities. As we assume in our physical model in section IV. B., stopping only 
takes place after the ball velocity has dropped below a certain value (in our case somewhere around 7 cm/s), at which 
point the rolling starts. This can be seen in Fig. |l^, where we have plotted the evolution of the x-velocity of the ball 
for the same cases as in Fig. The point where the ball starts to move in bounces much smaller than a ball radius 
(the stage prior to rolling) is marked with a small arrow for each trajectory. 

It can also be seen from these curves that after a very rapid drop in the velocity after the first few collisions with 
the plane, the velocity seems to decrease in a linear fashion, though with a slope that seems to depend slightly on the 
initial velocity. Then, when the particle enters the phase of the motion consisting of a number of bounces with each 
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ball it passes (and eventually some rolling), the friction force experienced by the particle drops to a much lower value 
and seems to be independent of the initial velocity. 

VI. CONCLUSION 

We have investigated by means of experiments and numerical simulations the movement of a single ball on a rough 
inclined plane, at small inclination angle. This has given us a better understanding of the energy dissipation at small 
angles. We have shown that a ball which has a large enough initial velocity Vq first bounces on the rough surface 
and suffers a constant friction force. Clearly in this case the ball cannot be trapped if its velocity is larger than the 
crossover velocity Vi. When the velocity reaches Vi, the friction force suddenly becommes viscous: the dynamics of the 
motion is now similar to that observed in regime B. The key for understanding these two mechanisms of dissipation, 
i.e. friction forces, is the difference in the nature of the collisions when the velocity is above or below Vi as explained 
above. We have also shown that the geometry of the surface plays an important role in the trapping of the ball. In 
regime A the ball is first slowed down gradually and when the velocity finally reaches a threshold value, Knm, (which 
appears to be independent of initial conditions) the ball is trapped. The trapping probability decreases linearly with 
the inclination of the plane. For the transition angle Ot, this trappping mechanism disappears and the ball crosses 
over into the dynamic regime B where it moves on the plane with a constant mean velocity and is subjected to a 
viscous friction force pjlCj]. In this regime, the ball can still be trapped by the occasional big defect on the surface 
but its trapping probability decreases exponentially with the angle of inclination. The fact that in both regimes A 
and B the friction force is viscous just before the ball gets trapped emphasizes the important fact that the difference 
between the two regimes A and B is not of dynamic origin, but due to two different trapping mechanisms. 

The "divergence" of the median, L*, of the stopping distance (defined as the length for which 50% of the balls are 
trapped) at 9t is seen to indicate clearly the transition between the A and B regimes. The values of L* depend on the 
initial velocity but the divergence always occurs at the same angle 0t- This shows that the location of the transition 
is independant of initial kinetic energy. Numerical simulations gave us "microscopic" details of the motion (like the 
time between collisions) which agreed with and confirmed the experimental measurements and our explanations. 
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Figure 1. Histogram of the stopping distance for D — 2mm (a) 9 — 3°5, (b) 9 — 5° , Xi 2mm 



Figure 2. Decay constant of the stopping distance distributions a as a function of the rough surface inclination angle 9, for 
D = 2mm. 
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Figure 3. Same as in figure ^ but obtained from the stochastic model, [12] 



Figure 4. Median distributions, L* , as a function of 6 for {o)D = 10.3mm, {o)D — 3mm, {t>)D = 1.6mm 



Figure 5. Variations of L* with 0, for D = 3mm, and two initial velocities: (o)Vi « 5 cm/s and {•)Vi ~ 45 cm/s. 

Figure 6. Distributions of the stopping distance for D = 10.3mm, 9 = 2°, and different initial square velocities V^^ oc Xi, 
from left to right Xi = 0.5, 4.5, 10.5, 16.5cm 



Figure 7. Average stopping distance L as a function of Vi^ oc Xi, for 6 = 2°, (A)D = 10.3mm and {o)D = 4.7mm 



Figure 8. Average stopping distance L as a function of Vi in arbitrary units, for = 2°, (A)D = 10.3mm, {o)D = 7.1mm, 
(>)D = 6.3mm, {•)D = 4.7mm, {<)D = 3mm and {o)D = 2mm. In filled lines the calculated variation for the corresponding D 
values. 



Figure 9. Average stopping distances for D = 5, 9 = 1° . 

Figure 10. Distance covered by the ball between sucessive collisions with the plane for three different starting velocities (a) 
Vi = 20cm/s, (b) Vi = 30cm/s, Vi = 40cm/s. The diameter of balls on the plane is 1 mm. 

Figure 11. Total distance covered by the ball as a function of time for the same cases as in Fig. ^ Full line corresponds to 
(a), dashed line to (b), dot-dashed line to (c). 



Figure 12. Velocity of the ball (same as in Fig. hoh. Full line corresponds to (a), dashed line to (b), dot-dashed line to (c). 
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